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The forward-backward asymmetry in e + e _ annihilation into a quark-antiquark pair is considered 
in the double-logarithmic approximation at energies much higher than the masses of the weak bosons. 
It is shown that after accounting to all orders for the exchange of virtual photons and W, Z -bosons 
one is lead to the following effect (asymmetry): quarks with positive electric charge (e.g. it, d) tend 
to move in the e + - direction whereas quarks with negative charges (e.g. d, u) tend to move in the 
e~ - direction. The value of the asymmetry grows with increasing energy when the produced quarks 
are within a cone with opening angle, in the cmf, Oo ~ 2Mz / \fs around the e + e~ -beam. Outside 
this cone, at 6o <C 8 <C 1, the asymmetry is inversely proportional to 9 . 

PACS numbers: 12.38.Cy 



I. INTRODUCTION 

The standard theoretical description of e + e~ annihilation into hadrons at high energies starts with the sub-process 
of e + e~ -annihilation into quarks and gluons, which is then studied with perturbative methods. It is usually considered 
as mediated by the exchange of all electrowcak (EW) bosons: e + e~ — > 7*, Z, W — > qq+gluons. One of the most well- 
known and successful predictions of the Standard Model is the forward-backward asymmetry, which has been studied 
for many years both theoretically and experimentally, particularly around the Z boson [Q. Next linear colliders will 
explore e + e~ annihilation at very high energies, probing further the Standard Model and eventually looking for New 
Physics. As it is well known, pure QED also gives rise to a forward-backward asymmetry even at low energies, albeit 
small, due to interference of one-photon and two-photon exchange diagrams. This effect persists at asymptotically high 
energies due to the multiphoton contributions in higher orders in a. Such multiphoton contribution in e + e~ — > u + ji~ 
was studied in Refs. 0] in the double-logarithmic (DL) approximation (DLA). The annihilation process in pi was 
considered in the following two kinematic regions: 

(i) Forward kinematics, when, in the center of mass frame (cmf), the outgoing /i + ( fj,~ ) goes in the direction of 
the initial e + (e~). 

(ii) Backward kinematics, when the outgoing /i + (fJ>~) goes in the e~ (e + ) -direction. 

These kinematical configurations refer to the case when the initial positive (or negative) electrical charges do not 
change the direction after the scattering, or they are affected by a major - almost backward - deviation. It was shown 
in (see also the review El) that at high energies the radiative DL corrections to the Born amplitudes are quite 
different for the forward and the backward kinematics. As a result, the cross section of the forward annihilation 
dominates over the backward one and therefore there is a charge forward-backward asymmetry: positive muons tend 
to go in the e + -beam direction and negative muons rather follow the direction of e~ . 

In the present paper we generalize these results, accounting also for DL contributions of multiple W,Z exchanges. 
We calculate below the forward and backward scattering amplitudes for e + e~ -annihilation into qq -pair at energies 
much higher than Mz and estimate the dependence of the asymmetry on the scattering angle and on the total 
energy of the process. We account for multiphoton exchanges as well as for multi- W, Z exchanges to all orders in 
the EW couplings in DLA. To this aim we construct and solve in DLA an infrared evolution equations (IREE) for 
the backward and forward scattering amplitudes. In doing so we follow the approach of Refs (j, Recently this 
approach was usedH for studying the double logarithmic asymptotics for lepton-antilepton backward scattering, for 
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scattering angles 9 ~ n in the EW theory. Furthermore, in Ref. [|| it is considered the case where the initial (final) 
state consists of a left (right) lepton and its antiparticle. However, because there is an essential difference in the EW 
theory in the description of the processes involving the left and the right fermions, in the present paper we generalize 
those results and consider the forward and backward annihilation with both left and right initial and final states. 
We also study these processes in a wider angular region 6 <C 1. Although we use the Feynman gauge through this 
paper, all our results are gauge invariant. The paper is organized as follows. Instead of calculating amplitudes of 
e + e~ — > quarks directly, we find more convenient to operate with SU(2) - invariant amplitudes of a more general 
process, the lepton-antilepton annihilation into quark-antiquark pair. Then in Sect. 2 we introduce such invariant 
amplitudes and show their relation to the forward and backward amplitudes for e + e _ -annihilation. The IREE for 
the invariant amplitudes are constructed in Sect. 3 and solved in Sect. 4. The solutions are obtained in terms of 
the Mcllin amplitudes corresponding to collinear kinematics. The IREE for the Mellin amplitudes are obtained and 
solved in Sect. 5. Then, we define and estimate the forward-backward asymmetry in Sect. 6. Finally, we discuss our 
results in Sect. 7. 



II. INVARIANT AMPLITUDES FOR LEPTON-ANTILEPTON ANNIHILATION INTO qq. 

We are going to account for the DL effects of exchanging the EW bosons to e + e~ -annihilation into quark-antiquark 
pairs of different flavours. When multiple W - exchange is taken into account, the flavour of the virtual intermediate 
fermion state is not fixed, though the initial and final states of the annihilation are well-defined. Because the EW 
theory organizes all fermions into doublets of the left particles and right singlets, this suggests that is more convenient 
to calculate first the scattering amplitude of a more general process, the annihilation of a lepton and its antiparticle 
into a quark - antiquark pair, and only after to specify the flavour of the initial and the final states. This turns to be 
easier because the effects of the violation of the initial SU(2)U(l) symmetry are in many respects neglected within the 
double-logarithmic accuracy. On the other hand that also means that the DLA can be applied safely only when the 
energy of the annihilation is much higher than Mz , My/ . At such energies the propagators of the SU (2) - gauge bosons 
W a (a — 1,2,3), D a b(k) ~ S a b/k 2 . The propagator of the U(l) -gauge boson B is 1/fc 2 in the same approximation. 
The SU(2) vertices of the W a interaction with the left fermions are gt a , where t a are the SU(2) generators and g 
is the coupling, whereas the vertex of the interaction of the field B with the left and the right fermions is g'Y/2, Y 
being the hypercharge and g' being the coupling. As in the most general process both the initial and final particles 
can be left and/or right, we consider all these cases separately. 

In this Sect, we consider the general case of the annihilation of the left lepton l k (p\) belonging to the doublet (v, e) 
and the antilepton li{p 2 ) from the charge conjugated doublet into the left quark q k (p^) belonging to the doublet 
(u, d) and the antiquark q~i'(p' 2 ) from the charge conjugated doublet. Therefore, the scattering amplitude M of this 
process is 



M = h(-p 2 ) q v (pi) Mg, l\vi) Q k ' (-P' 2 ) (1) 
where the matrix amplitude M^, has to be calculated perturbatevely. 

We consider the kinematics where, in the cmf, both particles of the produced pair move close to the lepton-antilepton 
beam. It corresponds to two kinematics: 

(i) forward kinematics when 

-t = -(p' 1 ~ Pl ) 2 «s=(p 1 +p 2 ) 2 tt-u = -(p' 2 -p 1 f , (2) 

(ii) backward kinematics when 

-u = -(p' 2 - Pl ) 2 « S ={pi+ P2 ) 2 ~-t = -{p' 1 -pi) 2 . (3) 

Then, replacing in (Q) the lepton-antilepton pair by e~ , e + and the quark-antiquark pair by fi~~ , /i + respectively, 
we see immediately that the electric charge almost does not change its direction in the forward kinematics (0) while 
it's reversed in the backward kinematics (|^) . Obviously that does not apply strictly for the annihilation into quarks 
because the electric charges of u -quarks and d -quarks are different in sign. Therefore t -kinematics is "forward" for 
the annihilation into a dd -pair and at the same moment it is "backward" for the annihilation into uu quarks. We 
will come back to this definition of backward and forward kinematics later, when we shall discuss the annihilation 
into u u and d d pairs, but until then we refer to (^) as t - kinematics, and (||) as u - kinematics. 
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In order to simplify the isospin structure, it is convenient to expand the matrix Mil, into a sum, each term 
corresponding to some irreducible representation of SU{2). In the t - kinematics (J2J) , the initial t - channel state is 
l k (Pi)Qi'(Pi)- Obviously, 



i*(pi)»'(pi) 



-8 k 5 b + (5 k S b , - -S k 5 b ) 



l a (pi)q b (p'i) 



b k s b a +2j2(tc)ktc) b a 



l a (pi)qM) 



(4) 



where the first term corresponds to the scalar and the second one - to the triplet representation of 577(2). Eq. (Q) 
suggests the the representation 

M^ = (P )f k/ M Q + (P 1 )f k ,M 1 (5) 
where Mo,i(s, t) are scalar functions and the singlet and triplet projection operators correspondingly are : 



(6) 



Similarly for the u - kinematics (||), where the initial u -channel state is l k (pi)q k (—p'2), irreducible SU(2) 
representations are obtained by symmetrization and antisymmetrization, 



Mg, = (P_) k \,M_(u, s) + (P + ) k \,M + (u, s) , 



with 



(P-) 



A-A-' 



°k°k' °k °k' 



(P- 



+ )kk' 



Using the projectors Pj , j = 0, 1, — , +, the invariant amplitudes can be easily obtained: 



(7) 



(8) 



3 (Pj)w'(Pj)k' k ' ' 

In the Born approximation, the amplitudes Mj defined in Eqs. can be written as 



(9) 



Mi 



R 



RAf orn {s) , Af orn (s) 



S + l€ 



"3 > 



«(-pa)[(l-7b)/2]7 v [(l+'7is)/2]«(pi)u(pi)[(l-7is)/2]7 v [(l+75)/2]«(-// 2 ) 



(10) 



where R denotes the normalized spinor factor and Aj orn are scalar functions of s, differing only in constant group 
factors aj . 

As we discuss the particular case of the left fermions we can drop the factors [(1 ± 7s)/2] and use the following 
definition: 



R 



u(-P2hMpi)u(p'ihvu(-p' 2 ) 



(11) 



For the left particles, the lepton and the quark hypercharges are Yi = — 1 and Y q = 1/3 respectively. The group 
factors aj for the Born amplitudes are 



a 



2g 2 +g' 2 Y l Y q 



a 1 



-gi + g'^Y, 



-3g 2 + g ,2 YiY q 



a+ 



g 2 + g' 2 YiY q 



(12) 
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The contributions proportional to Y\Y q in Eq. J12] ) come from the Born graph where the lepton line is connected 
to the quark one by the B -field, the other contributions come from the Born graphs with propagators of Wj -fields. 
Accounting for all DL corrections transforms the coefficients a,j into invariant amplitudes Mj, 

Mj = RMj(s,u,t) , (13) 

where, in DLA, Mj depend on s, t, u through logarithms. When Mj are calculated, Eqs. (^,||) allow us to 
immediately express the amplitudes e + e~ — * uu and e + e~ — > dd in terms of Mf. 



M(e + e~ -> m) = RMx(s,t), 

M(e + e- ^dd) = R [M a (s, t) + M x {s, t)] /2 , (14) 



for the annihilation in t -kinematics (0). Similarly, 



M(e+e" -» uu) = # [M_(s, t) + M+(s, <)] /2 , 

M{e + e~^dd) = RM + (s,t) , (15) 



for the annihilation in u -kinematics 



III. EVOLUTION EQUATIONS FOR THE INVARIANT AMPLITUDES M j 



In this section we calculate Mj in the high energy limit, by constructing and solving an IREE for them, as a 
generalization of the evolution equations derived earlier in QCD. This approach exploits the evolution of scattering 
amplitudes with respect to the infrared cut-off fi in the transverse momentum space. Transverse momenta of all 
virtual particles are supposed to obey 



ki± > n, h± ±px,P2- (16) 

The value of the cut-off \x must not be smaller than any of the involved masses, otherwise it is arbitrary. Introducing 
\x makes also possible to neglect the masses of all involved quarks and to restrict ourselves to consider the evolution 
of Mj with respect to [i only. Then one can take in the final formulas /i of order of the largest mass involved. In DLA 
we can also neglect the difference between the masses of EW bosons Mw and Mz, putting in the final expressions 



/i = M = M z w M w . (17) 

First we consider the annihilation in t -kinematics and construct the IREE for Mj with j = 0,1. According to (Q), 
t is small compared to u, —s. To bound it from below we assume that 



s > -t > (j, 2 . (18) 

The main idea of the IREE consists in evoluting the invariant amplitudes with respect to the infrared cut-off /x by 
applying to them the differential operator 



Vd/V , (19) 

in the form 

2 dM i dMj dMj 

-^V = ^ + *f' (20) 

where we use u w — s in this forward kinematics and have introduced the notations 



p = ln(s/ M 2 ), 7? = ln(-t/V 2 ) 



(21) 
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In order to obtain the rhs of last equation, we have to take into account the factorization of DL contributions of 
virtual particles with respect to p, where p is the lowest limit of integration over k±. In turn, this minimal k± acts 
as a new cut-off for other virtual momenta (see pi || for details). When the virtual particle with the minimal k± is 
a EW boson, one can factorize its DL contributions as shown in Fig. [I]. Applying then the differentiation (19) and 
the projection operators Pj we obtain with the help of Eq. (||) the contributions Go, G± to the EW singlet and triplet 
parts of the IREE respectively. 

Before writing the explicit expressions for Go.i, we want to discuss the general structure of these contributions. 
Integration over longitudinal momentum of the factorized boson with momentum k in graphs (a) and (b) in Fig, |l| yields 
ln(— s/kjj whereas the same integration in graphs (c) and (d) yields ln(— u/k^J w \n{s/k\). Similarly, graphs (e),(f) 
yield ln(— t/kj_). Besides these logarithms, each graph in Fig. [j] contains also an integration of Mj{s/k 2 \_,t/k 2 ± )/k 2 _ 
over k"j_, with the lowest limit p 2 . After differentiation (Il9|), we arrive at 



c,j 4 J 5,F).. k (?)-w..n? 

where the the quantities (b s ) , ho i will be explicitly given later 

V / o.i 



s t 

9 7 F 



Let us introduce 



->(±) 



In 



/'- 



±ln 



(23) 



so that ln(— s/p 2 ) = p^ + p( ■* and ln(— u/p 2 ) = — p^ \ Obviously, p^ and p( ■* are symmetrical and 
antisymmctrical functions with respect to replacing s by u. It is convenient also to introduce the invariant amplitudes 



M^{> with the same properties 



, l 



o,i 2 



so that Mq,i = Mqi + Mq x . Then for signature ampltudes G^ defined as 



(24) 



Gg = i[Go,i(M)±Go,i(u,t)] 



from Eq. (22) we obtain the following expressions 



(25) 



4? - g^E [(^C" (+, «! +, + (» ,J, )l>" M ! _ ' 

3=0,1 1 

<*> = f(" u, )r> (+,M ^ ,+ ( !>&i )i> HM i +i 



3=0,1 



^2^0,1??M ( J , 



(26) 



where 



( b (i)) (±) = (6«) T (ftW) . (27) 

V /0,1 V / 0,1 V /0,1 

Besides an EW boson, in kinematics (|) a i -channel virtual fermion pair, as shown in Fig. ||a, could also attain 
the minimal transverse momentum. However, DL contributions arising from the integration over this pair momentum 
could only come from the region k\ > —t 3> p 2 ■ Hence they do not depend on // in kinematics ([l8]) and must vanish 
when differentiated with respect to p. The same is true for the Born amplitudes (fio|). 



1 In the Regge theory, amplitudes are called the positive and negative signature amplitudes, and we use these notation below. 
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As soon as ln(— s/p 2 ) = ln(s/// 2 ) — iir, in kinematics (pj) with p = ln(s//i 2 ) w ln(— u/p 2 ) we obtain 

„(+) 



Z7T 

p - — sign(s) 



P l J = -ysign(s) 
It is assumed in DLA that \n(s/p 2 ) 3> 7r. This means that 



(28) 



(29) 



Similarly, in DLA in each order of the perturbative expansion the amplitudes Af- ' dominate over Mj } for one 
power of ln(s//i 2 ). By the same reason the amplitudes M- are mainly real, and we can assume 



Afj +) sa > 



r(+) 



M 



(-) 



(30) 



Combining Eqs. ( |2C|j2q) and using Eqs. (|3qj28|J 29|) leads us to the following IREE where negligible in DLA terms 
M^'p^' are dropped and terms M^'p^ ~ M^>p^' are retained: 









E K 

_J=0,1 


) (+) 


dp 


drj 


0,1 






= i p 


E K 

i=o,i 




dp 


dr] 


0,1 



r(+) 



\ 1 

~2/ 8^ 



E 

.7=0,1 



(31) 



C-) 



8tt- 



Let us proceed now to the u - kinematics (^). Using projection operators (|q) instead of (0) one can consider the 
annihilation in the it -kinematics and obtain the IREE for invariant signature amplitudes M_ , Afl introduced in 
a way similar to that one used for Mq\ '. As the amplitudes Afl , Afl correspond to SU(2) singlet and triplet 
representations, similarly to Mq\\ we can easily obtain IREE for u -kinematics from Eq. (31) with the replacement 



(32) 



Indeed, adding the restriction 



s > -u > /T 



(33) 



to Eq. (||), the derivation is quite similar to the previous one done for the t - kinematics thus leading to the same 
structure of the IREE for the amplitudes Afi^l. Therefore one can write down the same IREE for all invariant 
signature amplitudes M- , with j = 0, 1 ,"—","+" , generalizing Eq. (|3l| ) . 

The next significant simplification of ( |3l| ) comes from the explicit calculation of the group factors (fe^')^. It turns 
out that 



K +, -K +, -(^t'-(^) w ^ 



(34) 



and consequently the IREE for the positive signature amplitudes become linear homogenous partial differential 
equations. They can be written in a more simple general way: 



dM, 



(+) 



dM 



(+) 



dp 



1 

8k 2 



(35) 
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where bj , bj , and rf will be specified below, so that the only difference between the equations for different amplitudes 
comes from the numerical factors bj, hf. 



b_ = 



9'\Yi-Y q ) 2 
4 

g'\Y l + Y q f 



bl 



b + = 



+ g'\Y l -Y q f 



%g 2 + g l2 {Y l + Y q y 



(36) 



h 



h- = 



3 g 2 + g l2 YY q 
2 

3g 2 - g l2 YY q 



hi 



h+ = 



-g 2 + g'^Yg 
2 

-g 2 - g' 2 YY q 



(37) 



The IREE for the negative signature amplitudes Mj ' are also linear partial differential equations, but not being 



homogenous, they involve positive signature amplitudes through a non-zero matrix r 



33 ■ 



dM { r ] ^ dM J 



(-) 



dp drf 



1 

8^ 



— ITT 



E^" M 3> 



( + ) 



(38) 



In order to write down the IREE for all Afj ' and M^ 1 in the same way, we have used in Eqs. fl35|]38j ) the variable 
rf so that 



r(+) 



rf = T) = ln(— tj (j, 2 ) 



(39) 



for t - kinematics and 



rf = x = ln(~u/p 2 ) 



(40) 



for u - kinematics. 

The non-zero numerical factors Vjji in Eq. (|38| ) are 



l2 (Y + Y q ) 2 2 2 , a (Y - Y q ) 2 
roo =r u = g , r 01 = r_+ = Sg , r w = r + _ = g , r — = r ++ = g 



(41) 



IV. SOLUTIONS TO IREE FOR THE INVARIANT AMPLITUDES M 



(±) 



It is easy to check that a general solution to Eq. (pa) is 



M 



(+) 



<$>f\p-rf) exp[-^(C, 77')] 



(42) 



with 



1 

8^ 



bjCrj' + ibj+hj)^- 



(43) 



C = p — rf = const and is an arbitrary function. We can specify it by imposing the boundary condition 



(44) 
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where Aj + \p) are the amplitudes for the annihilation in the "collinear kinematics", i.e. in the kinematics where 
quarks are produced close to the direction of the beam of initial leptons, with the value of cither t or u much smaller 
than those fixed by Eqs. (p^ , |33| ). Of course one must use a separate boundary condition ( |iij ) for the t and u - 
kinematics. We define the notation "collinear t -kinematics" for 



and the notation "collinear u -kinematics" for 



-t < pf 



-u < p 



(45) 



(46) 



It is convenient to use the Mellin transform to represent signature amplitudes Ay 1 ' in the "collinear kinematics" 
IBpl): 



Af\ P ) 



(47) 



where 



(±) _ 



exp(— tirtu) ± 1 



(48) 



are the well known signature factors. At asymptotically high energy s the region of small uj, lu <C 1, is dominating 
in integral (|47|). This allows one to exploit the following approximations: 



t(+) 



VELO 



Eq. (|47]) implies that the positive signature amplitudes M- + ^ (p, 77') in the kinematic regions (|18|J33|) can be easily 



(49) 



expressed through the Mellin amplitudes Fj ' (to) : 

Ml + \p, V ')=exp[-cl> j (C, rj')} 



2m 



e("^')^(+'( u )f) +) ( w ) 



(50) 



On the other hand, as stated earlier, the IREE (BSfl for the negative signature amplitudes M- ^ are not homogeneous 



in contrast to Eq. (|35| ) Besides the amplitudes Afj \ they also involve the positive signature amplitudes Mj +> ■ In 
order to solve Eq. (Bq), we have to use again the boundary condition 



(+) 



?7'=0 



It is easy to check that the solution to the Eq. (38) satisfying Eq. (|5l| ) is 



(51) 



Mj ) (p,7 ? ')=exp[-^(C,ry')] 



4i } (p-r/)- / drexp[^(C,r)] 



-ITT 

~2~ J \ s- 



(52) 



Applying the Mellin transforms ( f47j ) and ( pp| ) for the amplitudes in Eq. ( |52] ) we can rewrite it through the Mellin 
amplitudes F- (w) as well. In order to simplify this procedure we use the substitution 
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which follows from the approximation (j49|) and is reasonable when the small u> region is dominating in the integral. 
Eventually we arrive at 



M<->(p,Tj') = exph^CT/)] 



2nt 



— y 



F ( +\lu) p 



33 



dr exp[<f>j (C, t) - <\> y (C, t)] 



(54) 



Eqs. ©H) show how one can obtain the amplitudes Mj when the Mellin amplitudes Fj (u>) are calculated. 



V. IREE FOR THE MELLIN AMPLITUDES F 



(±) 



In order to calculate -Fj ± ' ) we have to construct IREE for "collinear kinematics "amplitudes Ay 1 ' (ft). The IREE 
for them differ from the IREE for amplitudes M- considered in the previous section, by the following reasons. 



(i) The amplitudes Aj in the kinematics j45, 46) depend on p = ln(s//i 2 ) only, and some graphs in Fig. 1 do 
not yield DL contributions to IREE for them. In particular, graphs (e) and (f) with factorized t - channel 
virtual bosons do not contribute to IREE for A^ while graphs (c),(d) with factorized u - channel bosons do 
not contribute to IREE for The lhs of the IREE for A^ turns to -fjfdA^ jd[i 2 that in terms of the 
Mellin variable to results in wF^ (lu). 

(ii) The DL contributions of the graphs in Fig. || also depend on /j 2 and therefore do not vanish when differentiated 
with respect to \x. As these graphs are convolutions of two amplitudes, their contributions become simpler after 
applying the Mellin transform (|47]). Then the differentiation —(i 2 d/d[i 2 of these graphs leads to the following 



contribution to the IREE for the Mellin transforms -F„ 



_2_ 

8tt 2 



(55) 



with cj = 1 for j = 0, 1 and Cj = -1 for j = " - " , " 



(iii) Though at the first sight the Born amplitudes A^ orn of Eq. (10) do not depend on /i 2 it is necessary to replace 
them by ajs/(s — /i 2 + «e). This form explicitly tells that s cannot be smaller than [i 2 and also it makes the 
Mellin transform for the Born amplitudes to be correctly defined. The Mellin transforms for Born amplitudes 
are therefore aj/tu- Applying — /i 2 <9/9/! 2 to them results in multiplying by lu. Hence the contributions of the 



Born amplitudes to IREE are just the constant terms aj 



As a result we arrive to the following IREE for Mellin amplitudes F^ (lu) 



L0Fj +) (ij) 



8ir 2 



dijj 



loF\ \lo) = OJ- + -4 



bi 1 



d(wFj-\uj) 



8tt 2 oj 



du> 



^3 

8tt 2 



(56) 
(57) 



The coefficients a^bj, Cj and r^y are listed in Table | for t - kinematics and in Table |0] for u - kinematics. 
Solutions to Eq. (|5J) can be expressed in terms of the Parabolic cylinder functions D p : 



F {+) (to) - a J Dpj^MM 



Xj D P] (uo/\j) 



(58) 
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where 



8tt 2 



Pj 



(59) 



In contrast, solutions to Eq. ( |57| ) can be found only numerically. In QED the negative signature amplitudes for 

annihilation were solved in Ref. H. It is interesting to note that 6q = in the IREE 



for the forward amplitudes M^ 1 of e 



-annihilation, and the differential equations (|5^) for Me^' 



kinematics ( |45| ) turn into purely algebraic equations. This result was first obtained in Ref. g]. Later it was proved j| 
that the IREE for the (colourless) scalar components of the SU(3) negative signature amplitudes of quark-antiquark 
annihilation into another quark-antiquark pair are also algebraic and therefore can be easily solved. The processes 
mentioned above are the only known examples of solving IREE for negative signature amplitudes 2 . In all those cases 
the intercepts of the negative signature amplitudes are greater than those for the positive signature amplitudes, though 
the difference amounts only to a few percents. Equations for the negative signature amplitudes always involve the 
positive signature amplitudes. It has been observed^) in a QCD context, that these amplitudes can be approximated 
by their Born values with good accuracy. Such an approximation can help in solving Eqs. (|57|). We do not consider 
explicit solutions of Eqs. (pj) in the present paper. Instead, we consider below only contributions of amplitudes with 
the positive signature M- . Combining Eqs. (^0 58) and introducing variable x — w/Aj, we arrive at the expression 



M 



(+) 



P°° dx 

(p, T}') = a,j exp[-^ (p - rf, 77')] / — < 



D P Ax) 



(60) 



It is useful here to split <f)j defined by Eq. ( |43| ) and to combine its part depending on p — rj with the exponent of 
the integrand in Eq. (|6C|). Then changing the integration variable x to I, where x = I + AjT]' we finally obtain 



M { j + \p,'n') = a j exp 



b j + h j i_ 
8tt 2 2 



2vri 



DpAl + XjT]') 



(61) 



where, for the case of t - kinematics with j — 0, 1, rf = rj = ln(— t/p 2 ) and for the case of u - kinematics with 

7/ = x = M-V^ 2 ), 3 = ';-","+". _ _ 

The exponential factor in front of the integral in Eq. ( |6l| ) is of Sudakov type. Actually it is a product of the Sudakov 
form factors of the left lepton and of the left quark. As follows from Eqs. ( |36| , |37| ), it does not depend on j, i.e. is 
same for all invariant amplitudes, 



S = exp 




(62) 



It corresponds to DL contributions of soft virtual EW bosons and vanishes in the final expressions for the cross 
sections when bremsstrahlung of soft EW bosons are taken into account. Assuming this to be done we can omit such 
Sudakov factors. 

Until now we have discussed the annihilation 11 — *■ qq for the case when the both initial and final particles were 
left, i.e. the spinors in Eq. ( p"l| ) were actually [(1 + ys)/2]u. It is clear that applying the same reasoning it is easy to 
construct IREE for amplitudes in t and u -kinematics with right fermions. Solutions to such IREE can be presented 
in the same form of Eq. ( |6l| ) with j = RR for both right leptons and quarks, and j = LR when the initial leptons are 
left whereas the final quarks are right, or j = RL vice versa, with: 



g> YiY q \ I bRR l 



OiRR — 4 



Q>LR 



->LR 



g' 2 (Y lT Y q ) 2 



h RR = ± 
h L R 



9' YtY q 



_ 1 g' Y t Y q 
- 1 - 2 



(63) 



2 In the context of the EW theory, the IREE for the backward scattering amplitude with the negative signature was solved in 
irrealistic case of a complex value of the Weinberg angle. 



for the 



11 



where =F signs in b and ± signs in h correspond to t and u - kinematics respectively. It is worthwhile to remind here 
that in the numerical estimations when using Eq. (|(33|)) one has to substitute for Yi, Y q for right and left fermions the 
correct EW hypercharges: Y — 2Q for right fermions and Y = 2(Q — T3) for left fermions. The same formulae for the 
invariant amplitude M^ in the collinear u -kinematics, with \ — 0; can be also obtained from results of Ref. jfj. 



VI. FORWARD - BACKWARD ASYMMETRY 



In this section we consider in detail the e + e~ - annihilation into quark-antiquark pairs of different flavours in the 
configuration of forward and backward kinematics. According to the terminology introduced in Ref. 0], also discussed 
earlier in section |n|, the forward kinematics for e + (pi)e~ (P2) — > u (Pi)'^0 5 2) ~ annihilation corresponds actually to the 
u - kinematics defined by Eq. (0) because the electric charges of the electron and the u - quark are opposite. On the 
contrary, the forward kinematics for the annihilation of e + e~ —> dd corresponds to the t - kinematics, as we have 
defined in Eq. (^|). Similarly, the backward kinematics for e + (pi)e~ [p2) — > uip'-^uij)'^ -annihilation is a t - kinematics 
and that for e + e~ — > dd is the u - channel one. From the results we have obtained in section [H] it follows that the 
amplitude for forward e + e~ — ► uu -annihilation, is expressed in terms of amplitudes M_, M + of Eq. (pl|): 



whereas the backward amplitude for the same quarks is expressed through the amplitude Mi of Eq. (pll) : 



M* = Mi « m[ +) . (65) 



Similarly, the forward amplitude M? for e + e — > dd is 



while the backward amplitude for this process is 



Mf = M+ w i\4 +) . (67) 

By forward kinematics for e + e~ — » qq -annihilation we mean that quarks with positive electric charges, u and d, 
are produced around the initial e + - beam, in the cmf, within a cone with a small opening angle 6, 

1 > 9 > 9 - (68) 

By backward kinematics we means just the opposite - the electric charge scatters backwards in a cone with the 
same opening angles. 

The differential cross section daF for the forward annihilation is 

da F = da® [iMf | 2 + \M! I 2 ] = dff<°> [F u + F d ] , (69) 
and similarly, the differential cross section das for the backward annihilation is 



da B = dff(°) [|M„ B | 2 + \Mf | 2 ] = da^ [B u + B d ] , (70) 

where da^ stands for the Born cross section, though without couplings. It absorbs the factor |i?| 2 defined by 
Eq. (O) which cancels in the expression for the forward-backward asymmetry which we define as 
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where 



F = F U + F d , B = B U + B d . (72) 

Before presenting numerical results we would like to discuss the asymptotic behaviour of the asymmetry A. Let's 
discuss the angular region defined by (pq). The Regge factors (s/t) w , (s/u) u are equal to 9~ 2 " in this region. Eq. ( |6l| ) 
states that the asymptotics of all scattering amplitudes Mj + ^ is determined by the position of the rightmost zero of 
the Parabolic cylinder function D pj (z). The location of the zeros of D p , (z) depends on the value of Pj. In vicinity of 
the rightmost zero Za(j)j), D Pj {z) can be represented as 

D Pj {z) = D' p .(zo(pMz - zoipj)) (73) 
which, after substitution in Eq. (^l|) leads to the following expression for the high energy asymptotics of M ■ + ' : 



Mj +) cx 0- 2 ^Mp 3 ) (74) 

When pj < 0, all zeros of D pj (z) are in the left half of the complex z -plane and their real parts are negative, so 
that Mj (9) oscillates and decreases when 9 — ► 9q. We remind that all Xj are positive. 

When < pj < 1, the rightmost zero becomes real but still negative, so in this case again decreases when 

9 — > 6*o, though without oscillations. 

Finally, when pj > 1, the rightmost zero is real and positive, so that Mj (9) increases when 9 —> 9 . 
On the other hand, in the angular region of small 9 



9 > 9 (75) 

in Eq. (|74j) 1/9 2 must be replaced by s/fj 2 . Therefore, grows with s in the angular region (f75j) when pj > 1 

and decreases when pj < 1 . 

Basicaly, finding zeros of D p {z) functions is a rather tedious procedure. In order to simplify it, one can use the fact 
that, at integer p, 



D p ^(z)/D p (z) = V^Hp^z/^/Hpiz/V^) 



(76) 



where H p are the Hermite polynomials. Now, finding the zeros of a polynomial is easier that for D p {z). As Pj are 
expressed through a,j,bj,Cj (see Eq. (59)), with a,j, bj and Cj given in Table 1, it is easy to see that M^ d grow while 
M^ d fall when 9^0. Indeed, the small - 9 asymptotics of is controlled by the rightmost zeros of D p _ and D. 
Eq. (||) and Table || give p + < p_ and 



f r i 3.g- - g Y t Y q 

p u = max [p_, p+\ = p-= 2 w 25 . (77) 

9 [Yi + Y q y 



The small - 9 asymptotics of M d is controlled by the rightmost zeros of D po and D P1 . Eq. (|59|) and Table | give 
P\ < po and 



f r , 3, 9 2 - g'^Y, _ 

Pd = max [po, Pi] = Po= 2 « 6 . (78) 

9 (Yi + Y q ) 2 

Both u - quarks in Eq. (|n]) and d - quarks in Eq. ([78]) are left ones. Amplitudes M^, M d are controlled by D p 
with p^ = pi and p d = p+ respectively, and both p£ and p d are negative. Therefore, the asymmetry A, which 
is small at relatively large 9, grows up when 9 decreases down to 9q. When the produced quarks are within the 
"collincar" angular region ([75]), the asymmetry A does not depend on 9 but turns to depend on s. Besides, the 
forward amplitudes M^ d increase with s in this angular region whereas the backward amplitudes M^ d decrease. As 
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a result, the asymmetry A is s -dependent in the "collinear" angular region fl7q ) and increases up to 1 when s — > oo. 
Results of numerical calculations of A in "collinear kinematics" for different quark flavours are presented in Fig. ^. 
The curves below 100 GeV were calculated in pure QED with the infrared cut-off = 300 MeV for the case of light 
u, d and s- quarks and with the cut-off /i equal to quark mass for heavy quark flavours. Obviously, the value of 
the asymmetry is the same for the produced quarks and antiquarks (though they go in opposite directions). The 
energy region in vicinity of 100 GeV is close to the threshold of EW- boson production and therefore the DL radiative 
corrections at such energies come from virtual photon exchanges whereas radiative corrections involving the weak 
bosons are non-logarithmic. 

Experimentally, quarks manifest themselves through hadron jets. Then to observe the effect of the asymmetry one 
could measure asymmetry in production of a leading charged hadron Ti.^ (a meson M^ or a barion B^) inside a 
narrow cone along initial e~ (or e + ) beam: 



A{U) = da(H-) + da(H+) " (?9) 

As Mr = {UD} and B~ — {UUD} where U denotes any of u, c, t- quarks and D stands for any of d, s, b- quarks, 
one should expect a more rapid rise with the annihilation energy y/s for A(M) than for A(B). The reason is that at very 
high energies only U and D are produced in e~- beam direction whereas production of U and D- quarks is suppressed as 
the corresponding intercepts are negative. Both U and D- quarks can become constituent quarks of a leading Mr and 
cannot be constituent quarks of a leading M + . Non-constituent quark contributions due to fragmentation U — > M + 



and D — > M + are subtracted in the numerator of Eq. (79) and can be suppressed in the denominator by choosing the 



minimal momentum fraction x m i n in the definition of the leading hadron to be large enough. For the case of barion 
production a leading U can become a constituent quark of B~ , but now a leading D- quark can become a constituent 
of B + as well. So, a leading proton can also be observed in the direction of e~- beam. At asymptotically high energies, 
the contribution of U- quark will dominate over contribution of D- quark and A(B) will tend to unity, though slower 
than A(M) . However the numerical estimations for high but finite energies shown in Fig. |J demonstrate a sizable 
difference in the energy dependence between A(B) and A(M). It is worthwhile to note here that the phenomenological 
hadronization factors for mesons, W(U — > M~) = W(D — > M~) = W(U — > M + ) = W(D — > M + ) , are cancelled in 
Eq. ( fr9| ) for A(M) . But the barion asymmetry A(B) involves the ratio 



W{U^B+) 

R -W(D^B+) ■ m 

The barion curve in Fig. ^ corresponds to R — 2.5 which is in a good agreement with u and d- quarks fragmentation 
simulated with the JETSETjn) for x p > 0.3 and E quark = 200 2000 GeV. A more detailed study of hadron 
asymmetry with thorough account of nonleading and non-constituent contributions will be done elsewhere. 

The IREE ( |5r] , |57| ) also describe 11 -annihilation into another lepton-antilepton pair, with the lepton being left, 
providing Y q is replaced by the lepton hypercharge. Therefore, Eq. (|5^) can be used also for studying the forward- 
backward asymmetry for 11 — > VV -annihilation. Let us note that for such purely leptonic processes, the forward 
(backward) kinematics is identical to t (u) - kinematics. 



VII. DISCUSSION 



Forward-backward asymmetry in e + e _ annihilation into a quark-antiquark pair has been considered in the double- 
logarithmic approximation at energies much higher than the masses of the weak bosons, taking into account to all 
orders the exchange of virtual photons and W, Z -bosons. In deriving our results we have not considered other channels 
than the annihilation into a quark-antiquark pair. Of course, other channels could be present, which could change 
the cross section of the reaction, but will not contribute in DLA to the forward-backward asymmetry. Therefore we 
expect, at very high energies, that accounting for other channels will not affect the asymmetry. For example, besides 
EW effects, there are perturbative QCD corrections to both the forward and backward annihilation. In DLA, these 
corrections result into multiplying Mp and Mb by the same Sudakov form factor 



S, 



qcd = exp 



3tt 



hr 



(81) 



which suppresses the soft gluon emission. However, these Sudakov form factor effects cancel when the gluon 
bremsstrahlung is taken into account Q. The same is true for the soft photon emission. 
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In addition to the soft emission, there are also important QED and QCD corrections accounting for hard photon 
and gluon bremsstrahlung||. Indeed, the impact of this "harder" emission should be considered in more detail. Here 
we just notice that according to the results of Ref. B accounting for such emission changes the arguments of D p in 
Eq. ( |6l| ) but does not change the position of zq (p) and therefore it should not lead to big changes in our formulae for 
the forward-backward asymmetry. As stated earlier, the II —> VV -annihilation in the kinematics s 3> — u ~ M 2 and 
with chiralities of the fermions corresponding to the particular case j = LR of Eqs. (|6l] , |63| ) (corresponding to the EW 
singlet amplitude) was considered in Ref. M. The solution of the IREE for the positive signature amplitude of this 
process obtained in that analysis coincides with our corresponding case of Eqs. (61 6^). Also, for annihilation at low 
cm energies, when the asymmetry A in QED is due to the DL contributions from the multiphoton exchanges only, we 
agree with Ref.0. It is easy to obtain the expression for the forward and backward QED amplitudes from Eq. (|6l]), 



by putting g 



0, replacing g' by 47ra and substituting the hypercharges Yi, Y q by the electric charges Qi 
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TABLE I: The coefficients of IREE Eqs.(||,||) for i-kinematics. The angle 9 here is the Weinberg angle. 
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TABLE II: The coefficients for IREE Eqs.(B6|J5j) for ii-kinematics. The angle 9 here is the Weinberg angle. 
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FIG. 3: Charge asymmetry A for e + e~ — > for different quark flavors in the "collinear angular region". Curves 1-4 were 
obtained in pure QED (valid below threshold of EW bosons production): 

1 - for u- quark with the infrared cut-off fi = 300 MeV; 

2 - for d and s -quarks with the same cut-off; 

3 - for c- quark with the infrared cut-off /i = 1.4 GeV; 

4 - for b- quark with the infrared cut-off fj, = 4.5 GeV. 

Curves U and D correspond to any of u, c, t- quarks and d, s, b- quarks respectively and were calculated in DLA for EW theory 
with the infrared cut-off p = 100 GeV. 
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FIG. 4: Estimation of charge asymmetry A of leading charged hadrons in e + e annihilation: The curve M is for the meson 
asymmetry and the curve B is for the asymmetry of barions. 



